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CALIBRATION OF POSITIONING SYSTEMS 

This application claims priority under 35 U.S.C. §§ 119 and/or 365 to 
U.S. Provisional Application Serail No. 60/201,424 filed in the United States of 
America on May 3, 2000; the entire content of which is hereby incorporated by 
reference. 

FIELD OF INVENTION 

The present invention relates to the positioning of mobile stations within a 
cellular telecommunications network, and more particularly, to the calibration of 
positioning systems within a cellular telecommunications network. 

BACKGROUND 

Legal as well as commercial demands are driving a great interest in 
developing methods for positioning mobile stations (MS's). While mobile 
positioning may be achieved through the use of external systems, like Global 
Positioning Satellite (GPS) systems, mobile positioning can also be achieved by 
means internal to the cellular system. Several cellular positioning methods have 
been proposed, each of which relies upon measurements such as time of arrival 
(TO A), round trip delay, or angle of arrival of communication signals. Cellular 
positioning methods can be divided into uplink or downlink methods, i.e., whether 
the base station (BS) performs measurements on mobile station originating signals 
(uplink methods) or the mobile station performs measurements on base station 
originating signals (downlink methods). One example of a proposed downlink 
method is the Enhanced-Observed Time Difference (E-OTD) method. The E- 
OTD method is based on measuring the time-of-arrival (TO A) of bursts 
transmitted from a BS on its broadcast control channel (BCCH) carrier. 
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All location or positioning methods require knowledge of a number of 
parameters, some time- varying and others of more or less constant nature. These 
parameters can be TO A measurements, angle of arrival measurements, signal 
strength measurement or other parameters known to be used in mobile station 
positioning. For example, the E-OTD method requires that the location of the 
base stations and the real time difference (RTD) between transmissions from pairs 
of base stations are precisely known. Whenever these conditions are not fulfilled 
the location accuracy can be degraded. 

As discussed above, the Enhanced-Observed Time Difference is one 
example of a downlink positioning method. The E-OTD method is based on three 
quantities: the observed time difference (OTD); the real time difference (RTD); 
and the geometric time difference (GTD). These three quantities are related by 
the basic E-OTD equation: 



where OTD is the time difference between the reception of signals from two base 
stations, RTD is the timing or synchronization difference between two base 
stations, and GTD is the difference of propagation path length between two base 
stations. FIG. 1 depicts these quantities which are used in the E-OTD method. 

Referring to FIG. 1, if bursts from base station BS1 (e.g. serving base 
station) are received by a mobile station MSI at time tl and bursts from a 
neighbor base station BS2 are received by the mobile station at time t2, then the 
observed time difference is calculated as 



Likewise, if BS1 transmits its bursts at a time t3 and BS2 transmits its bursts at time 
t4, then the real time difference is calculated as 



OTD = RTD + GTD 



(1) 



OTD = t2 - tl. 



(2) 



RTD = t4 -t3. 



(3) 
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Finally, if dl is the propagation path length from BS1 to MSI, and d2 is the 
propagation path length between BS2 and MSI, then the geometric time difference 
is calculated as 



where dl and d2 are the distance of the propagation path from BS1 and BS2, 
respectively, and c is the speed of the radio waves, usually taken as the speed of 
light in vacuum. 

As illustrated in equation 5 below, if the OTD and RTD is known, the 
GTD can be determined in order to derive the location (i.e., position) information. 



Since the GTD is based on the difference in distance between a mobile station and 
two base stations, a known GTD defines a hyperbola upon which the Mobile 
station can be located. The position of the mobile station is given by the 
intersection of the hyperbolas defined by the GTDs, as shown in FIG. 2. 
Accordingly, at least two GTD values are required in order to locate the mobile 
station. 

Once a mobile station has measured the OTD values, in order to obtain the 
desired GTD values, the RTDs must be known. In addition, in order to obtain the 
location of the mobile station from the GTDs, the location of the base stations 
must be known. There are several methods for obtaining the RTD and the base 
station location. 

One method of obtaining the RTD values is to measure the real time 
differences using Location Measurement Units (LMUs). In its simplest form, an 
E-OTD LMU consists of an E-OTD capable mobile station placed at a fixed 
known location in the network. The E-OTD LMU measures the OTDs between 



GTD =(d2-dl)/c, 



(4) 



GTD = OTD - RTD 



(5) 
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pairs of base stations, and assuming that the locations of the LMUs and the base 
stations are known (i.e., the GTDs are known), RTD values can be calculated 
using the following equation: 

RTD = OTD - GTD (6) 

Depending on the network and the location of the LMUs (i.e., co-located at a base 
station site or at a location between base stations), the ratio between base stations 
and LMUs in the network is estimated to be 3:1 (i.e., for every 3 base stations in 
the network one LMU is required). Drawbacks of the LMU approach are that it 
requires new hardware to be deployed in the network, and depending on the 
stability of the base station clocks, the LMU's may have to report RTDs quite 
often, which increases the network load. 

An alternative method of determining the RTD values is to synchronize the 
network. This is the case, for example, in IS95 where GPS is used as a timing 
source. Also for GSM and TDM A, synchronization is being discussed and is 
indeed a requirement for deployment of a technology called compact EDGE. In a 
synchronized network, the base stations (or LMUs) don't need to report the RTDs 
since they are assumed to be constant. However, the synchronization is usually 
obtained in the digital baseband domain and not the air interface. It will be 
recognized that synchronization in the air interface is important for accurately 
determining location as a number of additional delay sources are present in the air 
interface which must be accounted for, including transmit filters, combiners and 
transmission lines. One solution for obtaining air interface synchronization is to 
use an LMU which measures the transmission on the air interface from its serving 
cell and relates the base station time to absolute time. The alternative of 
synchronizing the network partially overcomes the drawbacks associated with the 
LMU method, however; the synchronization needs to take the air interface into 
account. The required level of synchronization on the air is 5/xs for IS95, which 
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is enough for communication, but far too lax for location services (5/is 
corresponds to a traveled distance of 1500m). Factory calibration of all delay 
sources, including transmit filters, combiners and transmission lines is one 
possible solution to accounting for the air interface, but this may be prohibitive for 
cost reasons. 

The base station and LMU locations are normally not known precisely by 
the operator today, or are not accurate to the level required by location systems. 
This is due to the fact that for communication purposes exact locations are not that 
critical. To deploy location systems, therefore, it is necessary to measure 
precisely the BS and LMU locations, possibly by maintenance personnel equipped 
with differential GPS receivers. 

The problems with base station location calibration is that maintenance 
personnel need to visit each site and be equipped with differential GPS receivers. 
This may be costly and furthermore it may not be possible to obtain GPS coverage 
everywhere, e.g., inside buildings. In some regions in the world, differential GPS 
might not even be available. In these regions, the operator will have to rely on 
the accuracy provided by the standard service mode of GPS, which gives location 
accuracy of around 100m, 95% of the time. 

Accordingly, the invention provides a method for calibration of a 
positioning system within a telecommunications network which overcomes the 
drawbacks associated with known calibration methods. The method of the 
invention has the advantage that no new hardware is required, as it is entirely 
based on measurements which are made by multiple mobile stations. 

SUMMARY OF THE INVENTION 

As a solution to the above described problems, a method is disclosed that 
provides an estimation of the bias errors in the parameters that are relevant for 
mobile positioning. 
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According to embodiments of the invention, a method and system for 
determining mobile positioning is provided which accounts for bias errors in the 
relevant positioning parameters. This is achieved performing a plurality of 
measurements associated with a plurality of mobile stations and creating 
calibration parameters based on estimated mobile positions and the plurality of 
measurements. The estimated position of the mobiles is then updated using the 
estimated calibration parameters. 

According to another embodiment of the invention, a method for 
estimating bias errors in relevant positioning parameters is provided. 

BRIEF DESCRIPTION OF THE DRAWINGS 

These and other objects, features and advantages of the present invention 
will be readily apparent to one skilled in the art from the following written 
description, read in conjunction with the drawings, in which: 

Figure 1 illustrates the quantities which define OTD, RTD and GTD; 

Figure 2 illustrates the principles of the E-OTD location method; 

Figure 3 illustrates the steps of the bias estimation procedure; 

Figures 4A and 4B illustrate resulting location accuracy for different noise 
variances, using timing calibration in accordance with the present invention; and 

Figure 5 illustrates an exemplary embodiment of the present invention. 

DETAILED DESCRIPTION OF THE INVENTION 

In the following description, for purposes of explanation and not limitation, 
specific details are set forth, such as particular techniques, steps, circuits, circuit 
components etc. in order to provide a thorough understanding of the present 
invention. However, it will be apparent to one of ordinary skill in the art that the 
present invention may be practiced in other embodiments that depart from these 
specific details. In other instances, detailed descriptions of well-known methods, 
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devices, and circuits are omitted so as not to obscure the description of the present 
invention with unnecessary detail. 

These and other aspects of the invention will now be described in greater 
detail in connection with a number of exemplary embodiments. To facilitate an 
5 understanding of the invention, many aspects of the invention are described in 
terms of sequences of actions to be performed by elements of a communication 
system. It will be recognized that in each of the embodiments, the various actions 
could be performed by specialized circuits, by program instructions being 
ffi executed by one or more processors, or by a combination of both. Moreover, the 

]2 10 invention can additionally be considered to be embodied entirely within any form 

y\ of computer readable storage medium having stored therein an appropriate set of 

ijjj 

HI computer instructions that would cause a processor to carry out the techniques 

described herein. Thus, the various aspects of the invention may be embodied in 
many different forms, and all such forms are contemplated to be within the scope of 

GJ 15 the invention. 

O The calculation of bias error estimates according to an embodiment of the 

invention is formulated and solved in a generic framework, such that the 
procedure is applicable to various positioning methods. In the generic case, the 
measurement equation for determining the location of a mobile station can be 
20 defined as 

+ ** k = h.^N (7 ) 



where Y fc is a generic m-dimensional measurement related to the kth mobile 
station, b is an n-dimensional vector of unknown biases, z fc is the position of the 
kth mobile station represented by the column vector (x k y k ) T , and e k is an additive 
noise term with covariance P k . Equation (7) can be written in compact form as 

r=/(z,*)+« (8) 



f 
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with 

The goal is to estimate both the location z k and the bias vector b. This is achieved 
by first ensuring that z k and b can be uniquely estimated. 

Generally it is known that the number of equations must be larger than the 
5 number of unknowns. This is a necessary requirement, but as shown later, not 
sufficient. Accordingly, the present invention requires that mN > 2N+n-l. 
y3 For uniqueness, it is also required that 

"P fo v h J = Az 2 >*> 2 ) => z x = z a ,A 1 = b 2 . (9) 

UJ 

W\ As an example, consider the case of range measurements with uncertain locations. 



in 



To simplify notation, the coordinates are represented as complex numbers, i.e., 
10 z = x + iy . In this case , r becomes 

r ^z k fb^=abs(z k 'b^ (10) 

where abs() represents the elementwise absolute values of the quantity in brackets 
and b s denotes the unknown sensor locations. Now b s2 is selected such that 

ft^expC/Ox^+c) (11) 

where c is a complex arbitrary constant and O is an arbitrary rotation angle. This 
process amounts to a translation and rotation of the sensor positions. Then 

=abs(z-exp(-i<&yb +c) 

k s2 (12) 

=a&s(exp( -/c&XexpO'OXz^ +c) ~b ) 

=r(expOO)(z Jk +c),^ f2 ) 
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Hence the range measurement corresponding to the mobile station at z k and 
sensors at b s could equally well have been originated from the mobile station 
located at exp(i<D)(z k +c) and sensors located at b s2 . It can easily be verified that 
this holds for all measurements, i.e., f(z,b sl ) = f(exp(i<&)(z+c),b s2 ). Accordingly, 
the full bias vector can not be estimated. However, if, for example, one sensor 
location is accurately known and the orientation relative to one other sensor is 
known, then c and <& can be selected such that b s2 contains only n-3 unknown 
parameters. 

Now turning to the estimation of z u ...z N , b, this can be achieved in a least 
squares sense by minimizing the loss function defined in equation (13) using a 
Newton-Raphson approach. 

N l 

nz,b) = £ - (Y k -r(z k ,b)fp- k l (Y k -r(z k ,b)) (13) 

k-\ 

With 6=(z r ,...,z£,£ r ) r , one iteration step looks like 

e=e 0 -F'- 1 (e 0 )r , (e 0 ) (i4) 

In addition, with the first derivative of the loss function with respect to position 
defined as 

dV/dz k = -G T k P- k \Y t -r(z k ,b)) (15) 

where G k is the derivative of the function r(z,b) with respect to position, defined 
as 



G k =dr(z,b)/dz T 



(16) 
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and, with the first derivative of the loss function with respect to the bias vector 
defined as 

N 

dVldb=Y,H T k P-*{Y k -rtz k ,by) (17) 

k = l 

where H k is the derivative of the function r(z,b) with respect to bias, defined as 

H k =dr(z k ,b)/db T (18) 



The derivative vector of the loss function of equation (13) can therefore be written 
as shown in equation (19). 

''(0)=[(G^Vr(z^ (19) 
The second derivatives are 

dV/dz t dz T t =G r k P k 1 G k (20) 
dV/dz,dz^0 (21) 
dVldzdb T =G T P~*H , (22) 

N 



dVldbdb T =Y^ H T k P' k H k (23) 



t=i 



dVldbdz T =H T P~^G, (24) 

k k k k v ' 



The second derivative matrix can thus be written as 
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0 
0 



0 
0 



G T P l G G T P l H 

N M N iV N N 

lc - 1 



(25) 



-4 

hi 



Iteration steps: 

The second derivative matrix can be partitioned as 



A B 

D C 



(26) 



where 



O 
p 







0 


0 




GlP^H x 


v4 = 


0 




0 








0 


0 


G T P ~ X G 

N N N 




G l P N H N 



D = 



By using results on inverses of the partitioned matrices, it can be shown that 



A 


B 


-l 


U 


V 


D 


C 




w 


X 



(27) 



5 



Similarly, the first derivative matrix can be partitioned into 

V =[E T F T f 



(28) 
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with 

^=[(CfP 1 " 1 (F 1 -r(z 1 ,6)) r ,.,«?^P^(F JV -r(z JV ,6)) 7 '] r (29) 



and 

F-f:H T k P- k \Y k -r{z b ,b)) (30) 



This means that the update for b will be as shown in equation (31). 

b=b Q -WEXF (31) 

As shown in Figure 3, the procedure for estimating the bias error begins as step 
301 by determining z k assuming no bias, i.e, b=0. This is achieved by selecting 
z k to minimize V(z k ,b 0 ), where b 0 is the apriori value of the bias vector. Call 
these values z^. In this case the minimizing z w satisfies V'(z k0 ,b 0 )=0. Then, at 
step 303, an approximation of z k as a function of b is derived as discussed above. 
But since the vector E consists of V'(z ko ,b 0 ), k= 1 , . . . ,N it follows that E = 0, 
hence equation (31) reduces to that shown in equation (32). 

b=b 0 -XF (32) 

Then at step 305, the approximation derived in step 303 is solved for b. It is well 
known in the literature that an explicit expression for X can be defined as shown 
in equation (33), provided the inverses exist. 

X=(C DA l B) x (33) 
The final expression for equation (32) then becomes 



--K-(t< P k l H k -f:(HlP k l G k ^^ (34) 
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Finally, in step 307, z k is determined using the updated estimate of b derived in 
step 305. 

This procedure for estimating bias errors can be utilized in known or 
proposed positioning methods, for example, a E-OTD positioning method with 
timing calibration error. In this E-OTD positioning method the mobile station 
performs pseudorange measurements so that the basic measurement equation is 



with 



p(z)= ] /(x-xf +(y-y s ) 2 



(36) 



where c is the unknown timing calibration bias, d is related to the unknown transmit 
time, and v is a noise term with covariance N, and x s and y s of equation (36) are the 
base station coordinate vectors. The goal is to estimate c and z (and d) 
simultaneously. The unknown scalar d is usually eliminated by defining the 
differences of the pseudorange measurements, T, as 



1 -1 0 0 

10-10 • 

0 

10 0 -1 



(37) 



with the remaining parameters defined as b =Tc, r(z 9 b) = Tp(z) + b,c = Tv 9 
cov(e) =P = TNT T , then the resulting equation is exactly (7). Note that a 
consequence of the reformulation above is that b and not c will be possible to 
estimate. It is straightforward to show that in this special case, G of equation (16), 
becomes 
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G = TJ (38) 
where J is the derivative of p, defined as 





x-x 

n 


Pi 


P n 


y-y 1 


y-y„ 


P, 





(39) 



Furthermore, H of equation (18), in this case, the derivative of p with respect to b 
becomes simply the identity matrix I. 

H=I (40) 

To illustrate the above procedure, a small network consisting of 5 base 
stations was simulated. The base stations were arranged in a circle, and 300 mobile 
stations were randomly placed in the vicinity of the base stations. Pseudoranges p 
were estimated by computing the MS-BS distances, adding a random calibration 
bias which was the same for all mobile stations and an individual transmission time 
bias term. White noise was added to the pseudoranges p. The locations were 
estimated assuming zero bias and in the second step the calibration bias was 
estimated. The mobile station locations were then recalculated using the estimated 
bias terms. The resulting location accuracy is illustrated in Figures 4 A and 4B for 
different values of the noise variance and the magnitude of the bias errors. The 
results without timing calibration (curves 401 and 404), with timing bias elimination 
(curves 402 and 405) and, for comparison, from a bias free run (curves 403 and 406) 
are shown. As illustrated by curves 402 an 405, the bias is almost completely 
eliminated in both cases using the method in accordance with the invention. 
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This procedure can also be applied to E-OTD methods where there is 
uncertainty in the BTS locations. In this E-OTD positioning method (i.e., with 
BTS location uncertainty), the measurement equation is identical to the previous 
case except 

r(z k ,b) = Tp(z k ,b) (41) 

where 

P(z) = /(*"(* + b x )f +(y-(y s +b y )f (42) 



and b = bj>^ . As discussed above, it is not possible to estimate z k and b uniquely. 
To restrict the search space it is assumed that the location of the first two base 
stations are well known. In this case, the matrix G is also identical to equation (38), 
however; the matrix H is defined as 



H TK 



(43) 



where K is the derivative of p(z) with respect to b x and b y , defined as 



K = 



0 
0 



P 3 
0 



0 
0 

0 
0 

-(*-*) 



0 
0 

-(y-y 3 ) 



o 
o 



o 
o 

o 

0 

-(y-y n ) 



(44) 



Figure 5 illustrates the principle of operation of the invention in connection 
with timing calibration bias. With reference to Figure 5, the Location Services 
(LCS) application 503 sends a request for a location to a mobile location center, 
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which is illustrated in Figure 5 as position calculation function (PCF) 507. The PCF 
507 upon receiving the request sends an E-OTD measurement request message to a 
mobile station, illustrated in Figure 5 as E-OTD MS#k 509. The Mobile station 
performs the necessary E-OTD measurements (e.g., pseudoranges) and sends the 
results back to PCF 507 in an E-OTD Measurement Response message. The PCF 
507 computes the position of the Mobile station using, for example, equation (5). 
The resulting estimate is sent to the location application requesting the location. In 
parallel, the location estimate, along with the previous estimate of the calibration 
bias and mobile station reported pseudoranges, are sent with a calibration data 
message to the Timing calibration function 501. The Timing Calibration Function 
501 collects received location estimates from PCF 507 until sufficient data is 
available to perform calibration bias estimation according to equation (34). The 
resulting estimate is forwarded to PCF 507 in a Timing Calibration report. 
Periodically, or depending on the synchronization architecture, PCF 507 also 
receives RTD reports from Location Measurement Units (LMUs) or Base Stations 
(BS) 505. The PCF 507 uses initially the RTD reports in the position calculation, 
and later includes the updated calibration reports from the Timing calibration 
function 501. It will be recognized that Figure 5 represents functional elements 
used in accordance with exemplary embodiments of the present invention and that 
these functional elements may be contained in one node in a radio communication 
system, or may be distributed across several nodes. 

The invention has been described with reference to particular 
embodiments. However, it will be readily apparent to those skilled in the art that 
it is possible to embody the invention in specific forms other than those of the 
preferred embodiments described above. This may be done without departing 
from the spirit of the invention. 

Thus, the preferred embodiment is merely illustrative and should not be 
considered restrictive in any way. The scope of the invention is given by the 
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appended claims, rather than the preceding description, and all variations and 
equivalents which fall within the range of the claims are intended to be embraced 
therein. 



